A big data approach to evaluate receipt of optimal care in childhood cerebral palsy.

PURPOSE: Through automated electronic health record (EHR) data extraction and analysis, this project systematically quantified actual care delivery for children with cerebral palsy (CP) and evaluated alignment with current evidence-based recommendations. METHODS: Utilizing EHR data for over 8000 children with CP, we developed an approach to define and quantify receipt of optimal care, and pursued proof-of-concept with two children with unilateral CP, Gross Motor Function Classification System (GMFCS) Level II. Optimal care was codified as a cluster of four components including physical medicine and rehabilitation (PMR) care, spasticity management, physical therapy (PT), and occupational therapy (OT). A Receipt of Care Score (ROCS) quantified the degree of adherence to recommendations and was compared with the Pediatric Outcomes Data Collection Instrument (PODCI) and Pediatric Quality of Life Inventory (PEDS QL). RESULTS: The two children (12 year old female, 13 year old male) had nearly identical PMR and spasticity component scores while PT and OT scores were more divergent. Functional outcomes were higher for the child who had higher adjusted ROCS. CONCLUSIONS: ROCSs demonstrate variation in real-world care delivered over time and differentiate between components of care. ROCSs reflect overall function and quality of life. The ROCS methods developed are novel, robust, and scalable and will be tested in a larger sample.IMPLICATIONS FOR REHABILITATIONOptimal practice, with an emphasis on integrated multidisciplinary care, can be defined and quantified utilizing evidence-based recommendations.Receipt of optimal care for childhood cerebral palsy can be scored using existing electronic health record data.Big Data approaches can contribute to the understanding of current care and inform approaches for improved care.

Bering: joint cell segmentation and annotation for spatial transcriptomics with transferred graph embeddings.

Single-cell spatial transcriptomics such as in-situ hybridization or sequencing technologies can provide subcellular resolution that enables the identification of individual cell identities, locations, and a deep understanding of subcellular mechanisms. However, accurate segmentation and annotation that allows individual cell boundaries to be determined remains a major challenge that limits all the above and downstream insights. Current machine learning methods heavily rely on nuclei or cell body staining, resulting in the significant loss of both transcriptome depth and the limited ability to learn latent representations of spatial colocalization relationships. Here, we propose Bering, a graph deep learning model that leverages transcript colocalization relationships for joint noise-aware cell segmentation and molecular annotation in 2D and 3D spatial transcriptomics data. Graph embeddings for the cell annotation are transferred as a component of multi-modal input for cell segmentation, which is employed to enrich gene relationships throughout the process. To evaluate performance, we benchmarked Bering with state-of-the-art methods and observed significant improvement in cell segmentation accuracies and numbers of detected transcripts across various spatial technologies and tissues. To streamline segmentation processes, we constructed expansive pre-trained models, which yield high segmentation accuracy in new data through transfer learning and self-distillation, demonstrating the generalizability of Bering.

Guided construction of single cell reference for human and mouse lung.

Accurate cell type identification is a key and rate-limiting step in single-cell data analysis. Single-cell references with comprehensive cell types, reproducible and functionally validated cell identities, and common nomenclatures are much needed by the research community for automated cell type annotation, data integration, and data sharing. Here, we develop a computational pipeline utilizing the LungMAP CellCards as a dictionary to consolidate single-cell transcriptomic datasets of 104 human lungs and 17 mouse lung samples to construct LungMAP single-cell reference (CellRef) for both normal human and mouse lungs. CellRefs define 48 human and 40 mouse lung cell types catalogued from diverse anatomic locations and developmental time points. We demonstrate the accuracy and stability of LungMAP CellRefs and their utility for automated cell type annotation of both normal and diseased lungs using multiple independent methods and testing data. We develop user-friendly web interfaces for easy access and maximal utilization of the LungMAP CellRefs.

Schizophrenia-associated NRXN1 deletions induce developmental-timing- and cell-type-specific vulnerabilities in human brain organoids.

De novo mutations and copy number deletions in NRXN1 (2p16.3) pose a significant risk for schizophrenia (SCZ). It is unclear how NRXN1 deletions impact cortical development in a cell type-specific manner and disease background modulates these phenotypes. Here, we leveraged human pluripotent stem cell-derived forebrain organoid models carrying NRXN1 heterozygous deletions in isogenic and SCZ patient genetic backgrounds and conducted single-cell transcriptomic analysis over the course of brain organoid development from 3 weeks to 3.5 months. Intriguingly, while both deletions similarly impacted molecular pathways associated with ubiquitin-proteasome system, alternative splicing, and synaptic signaling in maturing glutamatergic and GABAergic neurons, SCZ-NRXN1 deletions specifically perturbed developmental trajectories of early neural progenitors and accumulated disease-specific transcriptomic signatures. Using calcium imaging, we found that both deletions led to long-lasting changes in spontaneous and synchronous neuronal networks, implicating synaptic dysfunction. Our study reveals developmental-timing- and cell-type-dependent actions of NRXN1 deletions in unique genetic contexts.

Spatial transcriptome profiling uncovers metabolic regulation of left-right patterning.

Left-right patterning disturbance can cause severe birth defects, but it remains least understood of the three body axes. We uncovered an unexpected role for metabolic regulation in left-right patterning. Analysis of the first spatial transcriptome profile of left-right patterning revealed global activation of glycolysis, accompanied by right-sided expression of Bmp7 and genes regulating insulin growth factor signaling. Cardiomyocyte differentiation was left-biased, which may underlie the specification of heart looping orientation. This is consistent with known Bmp7 stimulation of glycolysis and glycolysis suppression of cardiomyocyte differentiation. Liver/lung laterality may be specified via similar metabolic regulation of endoderm differentiation. Myo1d , found to be left-sided, was shown to regulate gut looping in mice, zebrafish, and human. Together these findings indicate metabolic regulation of left-right patterning. This could underlie high incidence of heterotaxy-related birth defects in maternal diabetes, and the association of PFKP, allosteric enzyme regulating glycolysis, with heterotaxy. This transcriptome dataset will be invaluable for interrogating birth defects involving laterality disturbance.

Single cell transcriptomic analysis of HPV16-infected epithelium identifies a keratinocyte subpopulation implicated in cancer.

Persistent HPV16 infection is a major cause of the global cancer burden. The viral life cycle is dependent on the differentiation program of stratified squamous epithelium, but the landscape of keratinocyte subpopulations which support distinct phases of the viral life cycle has yet to be elucidated. Here, single cell RNA sequencing of HPV16 infected compared to uninfected organoids identifies twelve distinct keratinocyte populations, with a subset mapped to reconstruct their respective 3D geography in stratified squamous epithelium. Instead of conventional terminally differentiated cells, an HPV-reprogrammed keratinocyte subpopulation (HIDDEN cells) forms the surface compartment and requires overexpression of the ELF3/ESE-1 transcription factor. HIDDEN cells are detected throughout stages of human carcinogenesis including primary human cervical intraepithelial neoplasias and HPV positive head and neck cancers, and a possible role in promoting viral carcinogenesis is supported by TCGA analyses. Single cell transcriptome information on HPV-infected versus uninfected epithelium will enable broader studies of the role of individual keratinocyte subpopulations in tumor virus infection and cancer evolution.

Informative missingness: What can we learn from patterns in missing laboratory data in the electronic health record?

BACKGROUND: In electronic health records, patterns of missing laboratory test results could capture patients' course of disease as well as ​​reflect clinician's concerns or worries for possible conditions. These patterns are often understudied and overlooked. This study aims to identify informative patterns of missingness among laboratory data collected across 15 healthcare system sites in three countries for COVID-19 inpatients. METHODS: We collected and analyzed demographic, diagnosis, and laboratory data for 69,939 patients with positive COVID-19 PCR tests across three countries from 1 January 2020 through 30 September 2021. We analyzed missing laboratory measurements across sites, missingness stratification by demographic variables, temporal trends of missingness, correlations between labs based on missingness indicators over time, and clustering of groups of labs based on their missingness/ordering pattern. RESULTS: With these analyses, we identified mapping issues faced in seven out of 15 sites. We also identified nuances in data collection and variable definition for the various sites. Temporal trend analyses may support the use of laboratory test result missingness patterns in identifying severe COVID-19 patients. Lastly, using missingness patterns, we determined relationships between various labs that reflect clinical behaviors. CONCLUSION: In this work, we use computational approaches to relate missingness patterns to hospital treatment capacity and highlight the heterogeneity of looking at COVID-19 over time and at multiple sites, where there might be different phases, policies, etc. Changes in missingness could suggest a change in a patient's condition, and patterns of missingness among laboratory measurements could potentially identify clinical outcomes. This allows sites to consider missing data as informative to analyses and help researchers identify which sites are better poised to study particular questions.

Outpatient hospital utilization after single event multi-level surgery in children with cerebral palsy.

PURPOSE: This study aimed to examine outpatient hospital utilization (number of specialties seen and number of visits to each specialty) in the year after single event multi-level surgery (SEMLS) in children with cerebral palsy (CP), and to determine if utilization differs across the medical center in the year after compared to the year before SEMLS. METHODS: This retrospective cross-sectional study used electronic medical record data of outpatient hospital utilization in children with CP who underwent SEMLS. RESULTS: Thirty children with CP (Gross Motor Function Classification System Levels I-V, mean age of 9.9 years) were included. In the year after surgery, a significant difference (p = 0.001) was found for the number of specialties seen, with non-ambulatory children seeing more specialties than ambulatory children. No statistically significant difference was found between the number of outpatient visits to each specialty in the year after SEMLS. Compared to the year before SEMLS, fewer therapy visits occurred in the year after SEMLS (p < 0.001) but significantly more visits to orthopaedics (p = 0.001) and radiology (p = 0.001). CONCLUSION: Children with CP had fewer therapy visits but more orthopaedic and radiology visits the year after SEMLS. Nearly half of the children were non-ambulatory. Examination of care needs in children with CP undergoing SEMLS is justified with consideration of ambulatory status, surgical burden, and post-operative immobilization.

Gross motor function prediction using natural language processing in cerebral palsy.

AIM: To predict ambulatory status and Gross Motor Function Classification System (GMFCS) levels in patients with cerebral palsy (CP) by applying natural language processing (NLP) to electronic health record (EHR) clinical notes. METHOD: Individuals aged 8 to 26 years with a diagnosis of CP in the EHR between January 2009 and November 2020 (~12 years of data) were included in a cross-sectional retrospective cohort of 2483 patients. The cohort was divided into train-test and validation groups. Positive predictive value, sensitivity, specificity, and area under the receiver operating curve (AUC) were calculated for prediction of ambulatory status and GMFCS levels. RESULTS: The median age was 15 years (interquartile range 10-20 years) for the total cohort, with 56% being male and 75% White. The validation group resulted in 70% sensitivity, 88% specificity, 81% positive predictive value, and 0.89 AUC for predicting ambulatory status. NLP applied to the EHR differentiated between GMFCS levels I-II and III (15% sensitivity, 96% specificity, 46% positive predictive value, and 0.71 AUC); and IV and V (81% sensitivity, 51% specificity, 70% positive predictive value, and 0.75 AUC). INTERPRETATION: NLP applied to the EHR demonstrated excellent differentiation between ambulatory and non-ambulatory status, and good differentiation between GMFCS levels I-II and III, and IV and V. Clinical use of NLP may help to individualize functional characterization and management. WHAT THIS PAPER ADDS: Natural language processing (NLP) applied to the electronic health record (EHR) can predict ambulatory status in children with cerebral palsy (CP). NLP provides good prediction of Gross Motor Function Classification System level in children with CP using the EHR. NLP methods described could be integrated in an EHR system to provide real-time information.

Hospitalizations Associated With Mental Health Conditions Among Adolescents in the US and France During the COVID-19 Pandemic.

Importance: The COVID-19 pandemic has been associated with an increase in mental health diagnoses among adolescents, though the extent of the increase, particularly for severe cases requiring hospitalization, has not been well characterized. Large-scale federated informatics approaches provide the ability to efficiently and securely query health care data sets to assess and monitor hospitalization patterns for mental health conditions among adolescents. Objective: To estimate changes in the proportion of hospitalizations associated with mental health conditions among adolescents following onset of the COVID-19 pandemic. Design, Setting, and Participants: This retrospective, multisite cohort study of adolescents 11 to 17 years of age who were hospitalized with at least 1 mental health condition diagnosis between February 1, 2019, and April 30, 2021, used patient-level data from electronic health records of 8 children's hospitals in the US and France. Main Outcomes and Measures: Change in the monthly proportion of mental health condition-associated hospitalizations between the prepandemic (February 1, 2019, to March 31, 2020) and pandemic (April 1, 2020, to April 30, 2021) periods using interrupted time series analysis. Results: There were 9696 adolescents hospitalized with a mental health condition during the prepandemic period (5966 [61.5%] female) and 11 101 during the pandemic period (7603 [68.5%] female). The mean (SD) age in the prepandemic cohort was 14.6 (1.9) years and in the pandemic cohort, 14.7 (1.8) years. The most prevalent diagnoses during the pandemic were anxiety (6066 [57.4%]), depression (5065 [48.0%]), and suicidality or self-injury (4673 [44.2%]). There was an increase in the proportions of monthly hospitalizations during the pandemic for anxiety (0.55%; 95% CI, 0.26%-0.84%), depression (0.50%; 95% CI, 0.19%-0.79%), and suicidality or self-injury (0.38%; 95% CI, 0.08%-0.68%). There was an estimated 0.60% increase (95% CI, 0.31%-0.89%) overall in the monthly proportion of mental health-associated hospitalizations following onset of the pandemic compared with the prepandemic period. Conclusions and Relevance: In this cohort study, onset of the COVID-19 pandemic was associated with increased hospitalizations with mental health diagnoses among adolescents. These findings support the need for greater resources within children's hospitals to care for adolescents with mental health conditions during the pandemic and beyond.

LungMAP Portal Ecosystem: Systems-Level Exploration of the Lung.

An improved understanding of the human lung necessitates advanced systems models informed by an ever-increasing repertoire of molecular omics, cellular, imaging, and pathological datasets. To centralize and standardize information across broad lung research efforts we expanded the LungMAP.net website into a new gateway portal. This portal connects a broad spectrum of research networks, bulk and single-cell multi-omics data and a diverse collection of image data that span mammalian lung development, and disease. The data are standardized across species and technologies using harmonized data and metadata models that leverage recent advances including those from the Human Cell Atlas, diverse ontologies, and the LungMAP CellCards initiative. To cultivate future discoveries, we have aggregated a diverse collection of single-cell atlases for multiple species (human, rhesus, mouse), to enable consistent queries across technologies, cohorts, age, disease, and drug treatment. These atlases are provided as independent and integrated queryable datasets, with an emphasis on dynamic visualization, figure generation, re-analysis, cell-type curation, and automated reference-based classification of user-provided single-cell genomics datasets (Azimuth). As this resource grows, we intend to increase the breadth of available interactive interfaces, supported data types, data portals and datasets from LungMAP and external research efforts.

CASK loss of function differentially regulates neuronal maturation and synaptic function in human induced cortical excitatory neurons.

Loss-of-function (LOF) mutations in CASK cause severe developmental phenotypes, including microcephaly with pontine and cerebellar hypoplasia, X-linked intellectual disability, and autism. Unraveling the pathological mechanisms of CASK-related disorders has been challenging owing to limited human cellular models to study the dynamic roles of this molecule during neuronal maturation and synapse development. Here, we investigate cell-autonomous functions of CASK in cortical excitatory induced neurons (iNs) generated from CASK knockout (KO) isogenic human embryonic stem cells (hESCs) using gene expression, morphometrics, and electrophysiology. While immature CASK KO iNs show robust neuronal outgrowth, mature CASK KO iNs display severe defects in synaptic transmission and synchronized network activity without compromising neuronal morphology and synapse numbers. In the developing human cortical excitatory neurons, CASK functions to promote both structural integrity and establishment of cortical excitatory neuronal networks. These results lay the foundation for future studies identifying suppressors of such phenotypes relevant to human patients.

CellDrift: inferring perturbation responses in temporally sampled single-cell data.

Cells and tissues respond to perturbations in multiple ways that can be sensitively reflected in the alterations of gene expression. Current approaches to finding and quantifying the effects of perturbations on cell-level responses over time disregard the temporal consistency of identifiable gene programs. To leverage the occurrence of these patterns for perturbation analyses, we developed CellDrift (https://github.com/KANG-BIOINFO/CellDrift), a generalized linear model-based functional data analysis method that is capable of identifying covarying temporal patterns of various cell types in response to perturbations. As compared to several other approaches, CellDrift demonstrated superior performance in the identification of temporally varied perturbation patterns and the ability to impute missing time points. We applied CellDrift to multiple longitudinal datasets, including COVID-19 disease progression and gastrointestinal tract development, and demonstrated its ability to identify specific gene programs associated with sequential biological processes, trajectories and outcomes.

International comparisons of laboratory values from the 4CE collaborative to predict COVID-19 mortality.

Given the growing number of prediction algorithms developed to predict COVID-19 mortality, we evaluated the transportability of a mortality prediction algorithm using a multi-national network of healthcare systems. We predicted COVID-19 mortality using baseline commonly measured laboratory values and standard demographic and clinical covariates across healthcare systems, countries, and continents. Specifically, we trained a Cox regression model with nine measured laboratory test values, standard demographics at admission, and comorbidity burden pre-admission. These models were compared at site, country, and continent level. Of the 39,969 hospitalized patients with COVID-19 (68.6% male), 5717 (14.3%) died. In the Cox model, age, albumin, AST, creatine, CRP, and white blood cell count are most predictive of mortality. The baseline covariates are more predictive of mortality during the early days of COVID-19 hospitalization. Models trained at healthcare systems with larger cohort size largely retain good transportability performance when porting to different sites. The combination of routine laboratory test values at admission along with basic demographic features can predict mortality in patients hospitalized with COVID-19. Importantly, this potentially deployable model differs from prior work by demonstrating not only consistent performance but also reliable transportability across healthcare systems in the US and Europe, highlighting the generalizability of this model and the overall approach.

International electronic health record-derived post-acute sequelae profiles of COVID-19 patients.

The risk profiles of post-acute sequelae of COVID-19 (PASC) have not been well characterized in multi-national settings with appropriate controls. We leveraged electronic health record (EHR) data from 277 international hospitals representing 414,602 patients with COVID-19, 2.3 million control patients without COVID-19 in the inpatient and outpatient settings, and over 221 million diagnosis codes to systematically identify new-onset conditions enriched among patients with COVID-19 during the post-acute period. Compared to inpatient controls, inpatient COVID-19 cases were at significant risk for angina pectoris (RR 1.30, 95% CI 1.09-1.55), heart failure (RR 1.22, 95% CI 1.10-1.35), cognitive dysfunctions (RR 1.18, 95% CI 1.07-1.31), and fatigue (RR 1.18, 95% CI 1.07-1.30). Relative to outpatient controls, outpatient COVID-19 cases were at risk for pulmonary embolism (RR 2.10, 95% CI 1.58-2.76), venous embolism (RR 1.34, 95% CI 1.17-1.54), atrial fibrillation (RR 1.30, 95% CI 1.13-1.50), type 2 diabetes (RR 1.26, 95% CI 1.16-1.36) and vitamin D deficiency (RR 1.19, 95% CI 1.09-1.30). Outpatient COVID-19 cases were also at risk for loss of smell and taste (RR 2.42, 95% CI 1.90-3.06), inflammatory neuropathy (RR 1.66, 95% CI 1.21-2.27), and cognitive dysfunction (RR 1.18, 95% CI 1.04-1.33). The incidence of post-acute cardiovascular and pulmonary conditions decreased across time among inpatient cases while the incidence of cardiovascular, digestive, and metabolic conditions increased among outpatient cases. Our study, based on a federated international network, systematically identified robust conditions associated with PASC compared to control groups, underscoring the multifaceted cardiovascular and neurological phenotype profiles of PASC.

Spaceflight Analogue Culture Enhances the Host-Pathogen Interaction Between Salmonella and a 3-D Biomimetic Intestinal Co-Culture Model.

Physical forces associated with spaceflight and spaceflight analogue culture regulate a wide range of physiological responses by both bacterial and mammalian cells that can impact infection. However, our mechanistic understanding of how these environments regulate host-pathogen interactions in humans is poorly understood. Using a spaceflight analogue low fluid shear culture system, we investigated the effect of Low Shear Modeled Microgravity (LSMMG) culture on the colonization of Salmonella Typhimurium in a 3-D biomimetic model of human colonic epithelium containing macrophages. RNA-seq profiling of stationary phase wild type and Δhfq mutant bacteria alone indicated that LSMMG culture induced global changes in gene expression in both strains and that the RNA binding protein Hfq played a significant role in regulating the transcriptional response of the pathogen to LSMMG culture. However, a core set of genes important for adhesion, invasion, and motility were commonly induced in both strains. LSMMG culture enhanced the colonization (adherence, invasion and intracellular survival) of Salmonella in this advanced model of intestinal epithelium using a mechanism that was independent of Hfq. Although S. Typhimurium Δhfq mutants are normally defective for invasion when grown as conventional shaking cultures, LSMMG conditions unexpectedly enabled high levels of colonization by an isogenic Δhfq mutant. In response to infection with either the wild type or mutant, host cells upregulated transcripts involved in inflammation, tissue remodeling, and wound healing during intracellular survival. Interestingly, infection by the Δhfq mutant led to fewer transcriptional differences between LSMMG- and control-infected host cells relative to infection with the wild type strain. This is the first study to investigate the effect of LSMMG culture on the interaction between S. Typhimurium and a 3-D model of human intestinal tissue. These findings advance our understanding of how physical forces can impact the early stages of human enteric salmonellosis.

Uncompensated mitochondrial oxidative stress underlies heart failure in an iPSC-derived model of congenital heart disease.

Hypoplastic left heart syndrome (HLHS) is a severe congenital heart disease with 30% mortality from heart failure (HF) in the first year of life, but the cause of early HF remains unknown. Induced pluripotent stem-cell-derived cardiomyocytes (iPSC-CM) from patients with HLHS showed that early HF is associated with increased apoptosis, mitochondrial respiration defects, and redox stress from abnormal mitochondrial permeability transition pore (mPTP) opening and failed antioxidant response. In contrast, iPSC-CM from patients without early HF showed normal respiration with elevated antioxidant response. Single-cell transcriptomics confirmed that early HF is associated with mitochondrial dysfunction accompanied with endoplasmic reticulum (ER) stress. These findings indicate that uncompensated oxidative stress underlies early HF in HLHS. Importantly, mitochondrial respiration defects, oxidative stress, and apoptosis were rescued by treatment with sildenafil to inhibit mPTP opening or TUDCA to suppress ER stress. Together these findings point to the potential use of patient iPSC-CM for modeling clinical heart failure and the development of therapeutics.

The genetic architecture of pediatric cardiomyopathy.

To understand the genetic contribution to primary pediatric cardiomyopathy, we performed exome sequencing in a large cohort of 528 children with cardiomyopathy. Using clinical interpretation guidelines and targeting genes implicated in cardiomyopathy, we identified a genetic cause in 32% of affected individuals. Cardiomyopathy sub-phenotypes differed by ancestry, age at diagnosis, and family history. Infants < 1 year were less likely to have a molecular diagnosis (p < 0.001). Using a discovery set of 1,703 candidate genes and informatic tools, we identified rare and damaging variants in 56% of affected individuals. We see an excess burden of damaging variants in affected individuals as compared to two independent control sets, 1000 Genomes Project (p < 0.001) and SPARK parental controls (p < 1 × 10-16). Cardiomyopathy variant burden remained enriched when stratified by ancestry, variant type, and sub-phenotype, emphasizing the importance of understanding the contribution of these factors to genetic architecture. Enrichment in this discovery candidate gene set suggests multigenic mechanisms underlie sub-phenotype-specific causes and presentations of cardiomyopathy. These results identify important information about the genetic architecture of pediatric cardiomyopathy and support recommendations for clinical genetic testing in children while illustrating differences in genetic architecture by age, ancestry, and sub-phenotype and providing rationale for larger studies to investigate multigenic contributions.

Betacoronavirus-specific alternate splicing.

Viruses can subvert a number of cellular processes including splicing in order to block innate antiviral responses, and many viruses interact with cellular splicing machinery. SARS-CoV-2 infection was shown to suppress global mRNA splicing, and at least 10 SARS-CoV-2 proteins bind specifically to one or more human RNAs. Here, we investigate 17 published experimental and clinical datasets related to SARS-CoV-2 infection, datasets from the betacoronaviruses SARS-CoV and MERS, as well as Streptococcus pneumonia, HCV, Zika virus, Dengue virus, influenza H3N2, and RSV. We show that genes showing differential alternative splicing in SARS-CoV-2 have a similar functional profile to those of SARS-CoV and MERS and affect a diverse set of genes and biological functions, including many closely related to virus biology. Additionally, the differentially spliced transcripts of cells infected by coronaviruses were more likely to undergo intron-retention, contain a pseudouridine modification, and have a smaller number of exons as compared with differentially spliced transcripts in the control groups. Viral load in clinical COVID-19 samples was correlated with isoform distribution of differentially spliced genes. A significantly higher number of ribosomal genes are affected by differential alternative splicing and gene expression in betacoronavirus samples, and the betacoronavirus differentially spliced genes are depleted for binding sites of RNA-binding proteins. Our results demonstrate characteristic patterns of differential splicing in cells infected by SARS-CoV-2, SARS-CoV, and MERS. The alternative splicing changes observed in betacoronaviruses infection potentially modify a broad range of cellular functions, via changes in the functions of the products of a diverse set of genes involved in different biological processes.